arXiv:1507.08382vl [hep-ph] 30M2015 


YITP-15-63, KUNS-2571 


Inhomogeneous chiral condensates and non-analyticity 
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We investigate inhomogeneous chiral condensates, such as the so-called dual chiral density wave of 
dense quark matter, under an external magnetic field at finite real and imaginary chemical potentials. 
In a model-independent manner, we find that analytic continuation from imaginary to real chemical 
potential is not possible due to the singularity induced by inhomogeneous chiral condensates at 
zero chemical potential. From the discussion on the non-analyticity and methods used in lattice 
QCD simulations, e.g., Taylor expansion, and the analytic continuation with an imaginary chemical 
potential, it turns out that information on an inhomogeneous chiral condensed phase is missed in the 
lattice simulations at finite baryon chemical potentials unless the non-analyticity at zero chemical 
potential is correctly considered. We also discuss an exceptional case without such non-analyticity 
at zero chemical potential. 

PACS numbers: 11.30.Rd, 21.65.Qr, 25.75.Nq 


I. INTRODUCTION 

Exploring the phase diagram of Quantum Chromody¬ 
namics (QCD) is one of the important and interesting 
subjects in nuclear physics, elementary particle physics, 
and astrophysics; for a review on the QCD phase diagram 
see Ref. [1], Ab initio lattice QCD simulations are a pow¬ 
erful and reliable method for analyzing the nonperturba- 
tive nature of QCD, but the well-known sign problem 
arises at finite real chemical potential (pr). To circum¬ 
vent the sign problem, several methods have been pro¬ 
posed so far: the Taylor expansion method, the reweight¬ 
ing method, the canonical approach, and the analytic 
continuation method from imaginary chemical potentials 
(/ii); see for example Ref. [2]. These methods are, how¬ 
ever, limited in the pr/T < 1 region, where T is the 
temperature. On the other hand, effective model calcula¬ 
tions are useful to investigate the phase diagram at finite 
T and pr, but it is still difficult to quantitatively discuss 
the phase structure of QCD due to a large ambiguity of 
the model. 

Inhomogeneous chiral condensed phases have been vig¬ 
orously studied in the QCD phase diagram at finite den¬ 
sity within chiral effective models, and also supported 
from Dyson-Schwinger studies of dense QCD (for a re¬ 
cent review see Ref. [3]), and their stabilities against ther¬ 
mal fluctuations have also been discussed [4, 5]. These 
phases have the quark condensate taking a spatially in¬ 
homogeneous configuration. One example with a one- 
dimensional modulation includes the so-called dual chiral 
density wave (DCDW) [6], where the chiral order param¬ 
eter is spatially modulated with a finite wavenumber ( q ) 
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in the single direction. Several studies based on mean- 
field calculations predict that the inhomogeneous phase 
can be realized in sufficiently high pr region. However, 
this situation could be changed if we introduce an exter¬ 
nal magnetic field. 

Magnetic aspects of QCD have been attracted much at¬ 
tention in the physics of heavy-ion collision experiments 
and compact stars (for a review see Ref. [7]). Recently, 
it has been reported that quark matter with the DCDW 
condensate could exhibit spontaneous magnetization [8], 
which may be related to the origin of strong magnetic 
fields in compact stars. In the presence of the external 
magnetic field, the most important and interesting point 
is that the DCDW phase appears even at small values 
of pr except pr = 0 within the massless Nambu-Jona- 
Lasinio (NJL) model [9-11]. This naturally raises the 
following question: Can we investigate inhomogeneous 
chiral condensates by using dense lattice QCD simula¬ 
tions at finite pr when the DCDW phase appears at 
very small pr! Here we assume that the inhomogeneous 
phase can appear in finite and discretized systems. 

In this paper, we discuss spatially inhomogeneous chi¬ 
ral condensates and the associated non-analyticity. Then 
we show the possibility that information on the inhomo¬ 
geneous phases can be missed in the lattice QCD simu¬ 
lations with the Taylor expansion method and the ana¬ 
lytic continuation method from the pi region to the pr 
region. We call this problem the “information missing 
problem.'’ We find that the applicable range of pr in the 
dense lattice QCD simulations at finite pr is strongly re¬ 
stricted, and such numerical simulations cannot describe 
the correct system with inhomogeneous chiral condenses 
if we neglect the non-analytic properties at p 2 = 0. Some 
problems in the canonical and reweighting methods are 
also discussed. In addition, an exceptional case where 
there is no non-analyticity at p 2 = 0 is presented. Here 
note that the same problem would exist in zero magnetic 
field, because the DCDW phase can appear also in the 
case without the external magnetic field. 
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II. INHOMOGENEOUS CHIRAL 

CONDENSATES AT IMAGINARY CHEMICAL 
POTENTIAL 

In realistic systems, the grand canonical partition func¬ 
tion (Z) must be real. It is also valid in the presence 
of inhomogeneous chiral condensates under a constant 
magnetic field, since the ground states of inhomogeneous 
phases correspond to solutions of Hamiltonian with ^r, 
i.e., Z(T,HYi,B\q) £ R, where B is the strength of the 
external magnetic field. This relation is simply mani¬ 
fested in the NJL model. The same discussion in the 
Polyakov-loop-extended NJL model requires some exten¬ 
sions of the formalism, e.g., considering a complex path 
contour [12-14] due to the model sign problem. 

Here we only focus on the DCDW phase, since it is well 
known that the so-called real kink crystal (RKC) phase 
does not appear in the fi\ region at least in the Gross- 
Neveu model [15]. Also, the RKC phase shows no signifi¬ 
cant dependence on the external magnetic field and thus 
appears at sufficiently high /tr, while the DCDW phase 
is extended to small /ir region. Therefore, we do not here 
treat the RKC condensate, since we are interested in the 
inhomogeneous phase at small jUr. 

Assuming that the DCDW type inhomogeneous con¬ 
densation arises in the presence of the uniform mag¬ 
netic field, the massless NJL Lagrangian density with 
the mean-field approximation takes the form [9] 

£ = ^^ 7 ^D At + /i7°-TO + 7 5 7 3 |)V'-^ ) ( 1 ) 

where ip is the quark field, 7 ^ are Dirac matrices, /1 is the 
quark chemical potential, G is the coupling constant, the 
constituent quark mass m = —2GA with A the constant 
amplitude of the DCDW, and the covariant derivative 
Dfj, = d,j, — ieA^ with e the electric charge and A M the 
external magnetic field. Here A M is directed along the 2 
axis, V x A = B z. The Lagrangian density given above 
is the one-flavor case but can be easily extended to the 
two-flavor one by considering the 7 r°-modulation. Then 
the quark determinant in the partition function takes the 
form 

In Z ~ Trln ( i^D M + qy 0 — m + ’j 5 j 3 q/2) . (2) 

Here we can show the symmetry of the system model- 
independently, while the analysis described above is 
model dependent. To this end, we evaluate Z by a 
Ginzburg-Landau (GL) expansion. 

In general, g-ocld terms are not banned in the expan¬ 
sion of Z with respect to q under the constant magnetic 
field: 

Z = a 0 + otiq + 0i2q 2 + &3Q 3 + oq ? 4 + ■ ■ • , (3) 

where a n (n = 0 , 1 , 2 ,...) are expansion coefficients de¬ 
termined by effective chiral models like the NJL model. 
Here the g-odd terms can be induced by a scalar quantity 


(q-B) without violating the rotational invariance. The 
coefficients a-i n (a 2 n+i) should be an even (odd) func¬ 
tion of B and q, respectively, which can be understood 
by taking the trace of Dirac space in Eq. (2). This im¬ 
plies that q of the DCDW can take an imaginary value 
in the finite qi region so as to manifest the reality of Z , 
i.e., Z(T,m,B;q = iqi) £ R where qi £ R. This sit¬ 
uation can arise in the case where an anomalous quark 
number density is proportional to q [10], since the quark 
number density (n g ) should be pure imaginary in the /q 
region and q acts as n q . A similar situation can also 
be expected, e.g., in complex chemical potential system, 
where Z should be a complex value. 

By summarizing the above discussions, we can reach 
the following scenarios: 

Scenario A — The DCDW phase does not appear at 
finite /q. 

Scenario B — The DCDW phase appears at finite /q, 
but the partition function becomes complex. 

Scenario C — The DCDW phase appears at finite /q, 
and the partition function is real. 

Each scenario is illustrated in Fig. 1. At a point very 
close to /i 2 = 0, it is enough to regard Z as only the 
first few terms in Eq. (3), since we can consider that q 
becomes quite small if there is no first-order transition at 
q 2 = 0. In the region of q 2 > 0, the stationary condition 
for q, dZ/dq = 0, yields q = — a\/2a 2 ~ q. Therefore, 
q 2 is described as a linear function of q 2 , as shown in 
Fig. 1. Here one can model-independently consider that 
Z contains a±, since it is known in Ref. [10] that the GL 
coefficient ol\ evaluated within the NJL model contains 
the Wess-Zumino-Witten (WZW) term [16] derived from 
chiral anomaly. In Scenario B, the wave number is real 
at finite /q, q £ R, while in Scenario C it becomes pure 
imaginary, iq £ R. Although we need a model calcula¬ 
tion to see which of scenarios is realized, our purpose in 
this paper is to qualitatively understand inhomogeneous 
chiral condensates and their associated non-analyticities, 
so that we leave it for future work. In the following, we 
discuss details of each scenario without model assump¬ 
tions. 


A. Non-analyticity and analytic continuation 

We discuss the analytic continuation from the /q to 
the q r region on the q 2 axis. Here we introduce two 
holomorphic functions, F\ and F 2 , which describe regions 
C x = {/i 2 £ R | 0 < (/i/T) 2 < e} and C 2 = {/r 2 £ R I - 
e < (/z/T) 2 < 0}, where e is a positive infinitesimal value. 
If F\ and F 2 have the connected domain at q 2 = 0, we can 
use the analytic continuation exactly. In Scenarios A and 
B, there exists the singularity at q 2 = 0 and hence the 
analytic continuation is no longer possible due to F\ ^ F 2 
for inhomogeneous chiral condensates. In Scenario C, on 
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FIG. 1: Schematic figures with respect to g 2 as a function of 
/x 2 at very close to /x 2 = 0 for each scenario. The /x 2 = 0 point 
becomes singular, except the solid lined case in Scenario C. 


the other hand, there is a possibility that the analytic 
continuation still works. 


1. Scenarios A and B 

In Scenario A, the /x 2 = 0 line becomes singular, so 
that the analytic continuation is impossible. This situ¬ 
ation, however, becomes non-trivial when a Taylor ex¬ 
pansion in terms of /xr at fx 2 = 0 is considered. The 
Taylor series cannot reproduce the non-analyticity and 
hence can describe only one side of the solution. To 
choice which solution is realized (whether the homoge¬ 
neous or inhomogeneous solutions), we need a restriction 
to the Taylor series, since we cannot mathematically de¬ 
termine the Taylor series just on such a singular point. 
In Scenario B, the situation is completely the same as 
Scenario A, where the /x 2 = 0 line is the singular line 
due to Im Z / 0 at finite fx The Taylor series with the 
restriction, Im Z = 0, may lead to solutions continued to 


the fin region. However, such a series is not mathemat¬ 
ically well-defined. We note here that one can consider 
a first-order transition at /i 2 = 0, where q 2 has a finite 
positive value in the /x 2 —>• —0 limit, but our conclusion 
is not changed in this scenario. 


2. Scenario C 

In Scenario C, the wavenumber has a nonzero real 
value, qr 7 ^ 0, in the C\ region, while in the C 2 region 
it becomes pure imaginary, q\ ^ 0. Here Z has only the 
q-even terms at fx 2 = 0 , since a 2 n+i is the odd function 
of /x in the expansion (3). Now we can further consider 
two situations, (I) and (II): 

Scenario C(I) — Assuming that q vanishes in the 
fj 2 —>• +0 limit, a^T 1 ,/x) should be positive at 
H = 0. In the fx 2 —> —0 limit, on the other hand, 
the coefficient of g 2 should be —02 < 0 for q = iq\. 
Thus, the point, q\ = 0, cannot be a minimum of 
Z , and qi does not vanish in the qi 2 —> —0 limit. 
In this case, there is the first-order phase transi¬ 
tion at /x 2 = 0. However, there is no reason that 
the translational (rotational) symmetry is sponta¬ 
neously broken in the [X 2 —> — 0 limit as well. 

Scenario C(II) — In Scenario C(I) it is assumed that 
a 2 (T,/j, = 0) ^ 0, but «2 = 0 can be acceptable. In 
this case, it seems possible to perform the analytic 
continuation because q 2 in the /x 2 —¥ ±0 limit can 
be smooth at /x 2 = 0. This situation is quite similar 
to the quark number density, n q , in the system with 
B = 0, since n 2 is positive (negative) at finite /xr 
(/ii) and n 2 is smoothly connected at fx 2 = 0 when 
we take the /a 2 —> ±0 limit. 

In Scenario C(I), the analytic continuation is impossible 
due to the first-order transition at fx 2 = 0. Although the 
Taylor series constructed just on the singular point is not 
mathematically well-defined, each solution of inhomoge¬ 
neous condensates may be reproduced under a suitable 
choice of input configurations. However, this situation is 
not realistic, since there is no reason to consider that the 
spontaneous translational (rotational) symmetry break¬ 
ing occurs at /j 2 = 0. Thus, we can exclude this situation. 
In Scenario C(II), on the other hand, there has the pos¬ 
sibility that the analytic continuation works because q 2 
becomes smooth at /i 2 = 0 as same as n 2 in the system 
with vanishing B. Only this scenario has the well-defined 
analytic continuation process and the Taylor series con¬ 
structed at /x 2 = 0. 

If inhomogeneous chiral condensates appear on the 
QCD phase diagram, the singular line, fx 2 = 0, should be 
exist. Note that this singular line can appear at smaller 
/xr /T than the natural boundary of the Taylor series in 
the absence of the magnetic field or the existence of the 
nonzero current mass. To avoid such a singularity in 
the analytic continuation, we can consider the complex 
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chemical potential. However, lattice QCD simulations 
and effective model calculations are quite difficult, since 
the singular region should have nonzero Im Z. 


B. Non-analyticity and dense lattice QCD 
simulations 

Now we discuss the consequences of the non-analyticity 
for several methods to work around the sign problem used 
in the lattice QCD simulations at finite chemical poten¬ 
tials. 

The trivial case is the analytic continuation method 
from /q to [i r [17-21]. This is nothing but our discussion 
about the analytic continuation. Thus, it is not feasible. 

The reweighting method [22-25] can go beyond the 
non-analytic point in principle. However, the overlap 
problem is important in this method. The /q region can¬ 
not be used for the creation of the probability in the im¬ 
portant sampling procedure, since in Scenario A the over¬ 
lap which can induce < 7^0 should be vanished or very 
small, and there is the sign problem at finite yq in Sce¬ 
nario B. Also, the p 2 = 0 point may have a small overlap 
problem. If we know regions where the sign problem does 
not arise and inhomogeneous chiral condensates exist, the 
reweighting method should completely work, while we do 
not know such a convenient region at present. 

The canonical approach [26-30] is based on the fact 
that the grand canonical partition function with pr can 
be constructed from the one with /q, which is valid if the 
quark number density is a good quantum number. In this 
method, we do not directly use the analytic continuation, 
so that the non-analyticity does not affect this method in 
principle. However, it is quite difficult to pick up infor¬ 
mation on inhomogeneous chiral condensates, since some 
inhomogeneous properties should be hidden when we cal¬ 
culate the grand canonical partition function or some ob¬ 
servables via the canonical partition function. This fact 
should be consistent with the results of a Lee-Yang zero 
analysis in QCD [31-33], since we only see the distribu¬ 
tion of zeros and hence it is difficult to clarify the exis¬ 
tence of inhomogeneous phases even if we can find any 
phase transitions from the behavior of zeros. 

On the other hand, the Taylor expansion method 
[34, 35] is nontrivial. Since this method is deeply related 
to the analytic continuation method, the Taylor series 
plays a crucial role. The Taylor series in the Taylor ex¬ 
pansion method constructed at q 2 = 0 is the holomorphic 
function. Thus, the /iR and /q regions cannot be simulta¬ 
neously reproduced because there is the non-analyticity 
at q 2 = 0. Mathematically, such Taylor series are not 
well-defined, but one can describe the both-side of solu¬ 
tions by considering the restriction to make the Taylor 
series describing the desirable limit by hand in the nu¬ 
merical calculations. This procedure may be affected by 
the numerical error. 

From the above discussion, we find that since the 
fj? = 0 point becomes the singular point, extreme care 


should be taken to perform the lattice QCD simula¬ 
tions with the Taylor expansion method, the reweighting 
method, the canonical approach, and the analytic con¬ 
tinuation method in the presence of the magnetic field. 
If the effect of the non-analyticity at q 2 = 0 are not 
correctly treated, information about the inhomogeneous 
chiral condensates should be missed in the Taylor ex¬ 
pansion and analytic continuation methods. Such a non- 
analyticity can be taken into account by imposing the 
restriction mentioned above. Therefore, the information 
missing problem should be carefully considered in the 
dense lattice QCD simulations. 

The information missing problem can be avoided in the 
complex Langevin method [36, 37] and the Lefschetz- 
thimble path integral method [38-40], since the dense 
lattice QCD simulation can be directly performed at fi¬ 
nite n r. However, both methods are far from perfection: 
there is no guarantee that the complex Langevin method 
leads to the correct answer when logarithmic terms ex¬ 
ist in the action (see, e.g., Refs [41-43], and Ref. [44] 
from the viewpoint of the singular drift term). The 
Lefschetz-thimble path integral method can provide the 
correct answer in principle, but the actual lattice calcu¬ 
lation method is not completed, e.g., how to include the 
multi-thimble contributions. 

The effective model approaches, on the other hand, are 
free from the information missing problem. The lattice 
QCD simulations without the assumption of inhomoge¬ 
neous chiral condensates can be possible by the suitable 
choice of spatial boundary conditions. If the effective 
models are constrained by such lattice QCD data at fi¬ 
nite /q, one can obtain the reliable model for the system 
without inhomogeneous phases. This method is noth¬ 
ing but the imaginary chemical potential matching ap¬ 
proach [45, 46]. However, information about inhomoge¬ 
neous chiral condensed phases can be introduced to the 
effective models by assuming the shape of the solution. 
This is a standard procedure to include the inhomoge¬ 
neous condensates in the effective model approach. By 
unifying the substitution procedure and the imaginary 
chemical potential matching approach, one can investi¬ 
gate full QCD. 


III. SUMMARY 

In this paper, we have investigated the properties of 
inhomogeneous chiral condensates at finite imaginary 
chemical potential in the presence of an external mag¬ 
netic field. From the reality of the partition function and 
the symmetry argument, we have considered two possi¬ 
ble scenarios for the DCDW at finite /q. In both sce¬ 
narios, the ff 2 = 0 point becomes the singular, which is 
induced by the inhomogeneous condensates. This singu¬ 
larity disturbs the investigation of inhomogeneous chiral 
condensates, and hence information of such condensates 
sometimes is missed, i.e., the information missing prob¬ 
lem. We found that the information missing problem can 
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appear in the lattice QCD simulation with the Taylor ex¬ 
pansion and analytic continuation methods. Problems to 
observe the inhomogeneous chiral condensates by using 
the reweighting method and the canonical approach are 
also discussed. Moreover, we considered the exceptional 
case where q 1 2 3 4 5 6 7 8 9 10 11 12 becomes negative. 

In the analytic continuation, holomorphic functions, 
e.g., the Taylor series, play a crucial role. If the fi 2 = 0 
point becomes the connected domain where two holomor¬ 
phic functions prepared in the [i,r and qi regions coincide, 
the present information missing problem does not arise. 
However, there exists the singularity at q 2 =0, which 
is induced from the different form of chiral condensates 
or the imaginary part of the partition function. Thus, 
the analytic continuation from the q 2 < 0 region is not 
possible, since the line /i 2 = 0 forms the singular line. 
This fact leads to that the Taylor expansion method re¬ 
quire an extreme care if the Taylor series is constructed 
at fx 2 = 0, the singular point. To correctly construct the 
Taylor series which continues to the /j,r region, some re¬ 
strictions are required in the numerical code. It should be 
noted that we can still consider Scenario C(II) which has 
the well-defined analytic continuation. However, even if 
which scenarios are realized, we must check the behavior 
of q 2 at fj, 2 < 0 before investigating the fi 2 > 0 region to 
clarify the existence of non-analyticity. 

For finite current quark masses, the singular line can 
be sifted at finite /ir region, which leads to the possibil¬ 
ity that the singular line appears at the smaller [Xr/T 
than the natural boundary of the Taylor series. In this 
case, the information missing problem becomes more se¬ 
rious. Consequently, the applicable range of the Tay¬ 
lor expansion method and that of the analytic continua¬ 


tion method are strongly restricted. It should be noted 
that we cannot estimate the position of the singular line 
by only using the lattice QCD data, unlike the natural 
boundary which can be estimated from the convergence 
behavior of the series. Meanwhile, a nonzero current 
quark mass gives rise to the configuration change of the 
DCDW condensate [47]. 

Possible promising methods to overcome the informa¬ 
tion missing problem at large /xr would include the unifi¬ 
cation of the imaginary chemical potential matching ap¬ 
proach and the inputting of the solutions of inhomoge¬ 
neous chiral condensates. In this case, one can obtain re¬ 
liable effective models at finite /ir without the inhomoge¬ 
neous phases, but the information on the inhomogeneous 
condensations can be restored by inputting the solution. 
In this method, the lattice QCD data without inhomoge¬ 
neous chiral condensates at finite fii is required but can 
be obtained if the spatial boundary conditions are suit¬ 
ably imposed. Therefore, we do not need details of infor¬ 
mation about the inhomogeneous condensates at finite 
fj\. Also, it is interesting to combine this approach with 
a nonperturbative method, e.g., the functional renormal¬ 
ization group method or the Dyson-Schwinger formalism. 
These topics will be considered in future studies. 
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